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Abstract. A phenomenological theory of exciton condensation in conditions of inhomoge- 
neous excitation is proposed. The theory is applied to the study of the development of an 
exciton luminescence ring and the ring fragmentation at macroscopical distances from the 
central excitation spot in coupled quantum wells. The transition between the fragmented 
and the continuous ring is considered. With assumption of a defect in the structure, a pos- 
sibility of a localized island of the condensed phase in a fixed position is shown. Exciton 
density distribution is also analyzed in the case of two spatially separated spots of the laser 
excitation. 
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1 Introduction 

Indirect excitons in double quantum well (QW) structures consist of electrons and holes driven 
to the separate wells by an external electric field. They are long living particles due to the 
hampered electron-hole (e-h) recombination. The enhanced lifetime (more than 3 orders of 
magnitude longer than that of the direct excitons in single QWs) makes it possible to create 
large concentrations of excitons in order to study effects of the exciton-exciton interaction, an 
exciton condensation m^JISllSISlinilZI, in particular. 

In AlGaAs and InGaAs based structures, a non-trivial feature in the photoluminescence 
spectra of excitons was observed jH El IS El IHl 13 • A spot of the laser excitation was reported 
to be surrounded by a concentric bright ring separated from the laser spot by an annular dark 
intermediate region. The distance between the ring and the spot grew with increasing intensity 
of the pumping and reached up to hundreds n m exceeding by much the exciton diffusion path. 
In certain cases, an internal ring was observed nearby the laser spot. 

At low temperatures (about 2 K), the external ring fragmented into a structure with a strongly 
evaluated periodicity on a macroscopic scale |2l IH lEl- Fragments followed the external ring, 
even when the excitation spot moved on the sample area or when the ring radius varied with 
irradiation intensity. With temperature rising, the external ring disappeared. 

Photoluminescence intensity of indirect excitons strongly grew in certain fixed spatial re- 
gions For any allocation of the excitation spot and its intensity, "localized bright spots" 
were observed, only when they were inside the area restricted by the external ring. These spots 
washed out with temperature growth. 

Under the excitation of the crystal by two spatially separated laser spots, it was revealed that 
in the process of rapprochement, the rings became deformed and slightly opened in the mutual 
direction and finally formed a common oval- like ring 
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The mechanism of the formation of the ring was suggested in [^1 [7] . It was based on two 
assumptions: 

1) in dark the well is populated with a certain density of electrons; 

2) holes are captured by the well with a larger probability than electrons. 

As a result, the irradiated structure develops two differently charged spatial regions. The laser 
spot and the region around are charged positively while the region far away remains negative. 
As electrons and holes can recombine only where they meet, the sharp luminescence ring forms 
on the boundary between the regions of the opposite charges. 

Assuming an attractive interaction between excitons, in the framework of a statistical app- 
roach papers jS[ IHI suggested an explanation of the fragmentation of the ring by means of 
the nucleation of spots of a condensed phase. Models of phase transitions in a system of indirect 
excitons were studied in ^01 • Finite lifetime of the excitons leads to the limitations on the size 
of the condensed phase areas, similarly to e-h droplets in bulk semiconductors. In a 2D system, 
these areas assume a shape of disk- like islands. The islands appear in the places of the highest 
exciton generation, i.e. on the ring. Due to the interaction between islands via the exciton 
concentration fields, there is a correlation in their positions eventually resulting in a periodicity. 

The statistical approach suggested in [SI [H] was based on the following assumptions: 

1) the exciton generation rate on the ring is a smooth function of coordinates, and is almost 
constant within the region of a condensed phase island; 

2) the distance between islands is much greater than the size of an individual island; 

3) the boundary between the gas and the condensed phase is sharp. 

In the present work we go beyond the stochastic approach of [HI E] and remove the above 
assumptions and solve the problem phenomenologically. We have incorporated a term containing 
free energy of excitons into the equation controlling the exciton density. The resulting equation 
was solved with account of the pumping and finite lifetime of excitons. This approximation 
allows treatment of inhomogeneous systems and obtaining new results, for instance, description 
of transition from a fragmented to a continuous luminescence ring with changing parameters of 
the system, development of an internal ring in additions to the internal one, an explanation of 
localized spots in the emission, a deviation of the shape of fragments from the spherical, the 
excitation of the system by two laser spots. Study of these essentially inhomogeneous problems 
in the framework of the statistical model would be very complicated. 

A mechanism leading to instability in exciton subsystem was suggested in the paper yet 
the problem of the structures of the exciton density forms outside the laser spot at a non-uniform 
pumping was not considered. 



2 The exciton generation rate on the ring 

In order to determine the exciton generation rate, we use the model of the ring formation 
suggested in jHlEj. We shall obtain solutions for the cases, which will be applied in the further 
study of the fragmentation. The equations for the electron and hole densities in a QW plane 
outside the laser spot are 

^ = D^An, + - Wn^n, - ^^1^, (1) 

dt Te 

^ = DhAn^, + Kh- WneUh - — . (2) 
dt Th 

Ue {nh) is the electron (hole) concentration, {Dh) is the electron (hole) diffusion coefficient, 
Te {Th) characterizes the time of the establishment of the electron (hole) equilibrium between the 
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Figure 1. The radial profile of the exciton generation rate G{p) in a QW plane. Pumpings are Ce = 
1.5 • 10^^ cni^^ s~^. = 3 • 10^^ cm"^ for the curves 1, 3, 4 and 1.33 times larger for the curve 2. 
W is 3 cm^ s^^ for the curves 1 and 2; 22 cm^ s~^ for the curve 3; 0.8 cm^ s~^ for the curve 4. 

well and the crystal outside, W is the e-h recombination rate, rig is the electron concentration 
in the well in absence of the irradiation. (Kh) is the electron (hole) creation rate, the radial 
profile of which is approximated by a Gaussian curve 



where q is the share of recombined electrons and holes, which participate in exciton formation 



The system of the equations and Q was solved numerically on a disk or on a rectangular 
plate with the sizes considerably exceeding the ring radius. The task was solved with the 
following initial conditions: ne{r) = n^, nh{r) = at t = 0, and the boundary conditions: 
rie (t) = rig, n/j (t) = at r ^ oo. We chose the following parameters of the system: Df, = 
200cm2s-\ Dh = 50cm2s-\ Tg = r/, = 10"^ s, n° = 10"^° cm'^, / =60^ m, q = 0.9. The 
values of D^, D^, Tg, coincide with the values in [Sj. 

The Fig. 1 shows the steady-state solutions of the equations and ((2)). Radial distribution of 
the exciton density has a sharp maximum at a certain distance from the center. This maximum 
moves farther away from the laser spot with increasing the pumping (see curves 1, 2 in Fig. 1). 

Width of the peak depends on the e-h recombination rate W. This rate is determined by 
the tunneling processes and depends on the width of the barrier which separates the QWs 
in the double well structure. With increasing W, the peak does not change its position, but 
narrows down and grows in height (compare curves 1, 3, 4 in Fig. 1). Its integral value remains 
approximately, within several percent, constant. 

3 Model of the system. Exciton density equation 

The microscopic theory of exciton condensation at inhomogeneous excitation is a difficult task, 
because we must take into account relaxation and finite exciton lifetime. Therefore, to build a 
solution, we apply a phenomenological approach. Our treatment of the problem is based on the 
following assumption: we assume that time of the post-excitation establishment of a quasi-local 
equilibrium of electron and hole densities and their binding into excitons is short, much smaller 
than time of the equilibration between different areas. The latter is related to slower diffusion 
processes. In this case, free energy of the quasi-equilibrium state can be considered as a function 




(3) 



The exciton generation rate is proportional to the product of rig and Uh- 



G = qWueUh 



(4) 



((/<!). 
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of the exciton density, which depends on the spatial coordinates. Obviously, such description is 
possible at macroscopic distances, much larger than the exciton free path, on which the wave 
function loses coherence. Then, a phenomenological equation for the exciton density can be 
written down as 

_-_ = -divj+G , (5) 

Ot Tex 

where n^x is the exciton density, Tcx is the exciton lifetime, j is the density of the exciton current: 

3 = -MVfi, 
M is the exciton mobility: 



M 



H is the chemical potential: 
6F 



6 Hex 

We chose free energy in the form suggested by the Landau model: 



K 2 

— {Vucx) + / (ricx) 



(6) 



Here, the term K (Vnex)^ /2 characterizes the energy of an inhomogeneity. Free energy density / 
is 

/ ("-ex) = ^ (ncx - n^cxf + I (nex - n^^xf + ^ (nex - <x)^ , (7) 

a, b, c, n^^ are phenomenological parameters, which can be obtained from the quantum mechan- 
ical calculations, or extracted from the comparison of the theory with experiment. 
After substituting equations © and ((TJ into equation the latter becomes 

^ = -A^nex + (3n^,2 _ ^^^c^ _ ^^^^ _ !^ + ^ - (Sn^ - /3) A^iex' + Anex^ (8) 



/3 



Here we introduce dimensionless units 

for length, time and exciton concentration, respectively. We assume that equation (jHj) holds also 
if the condensed phase is an e-h liquid. Then F is a function of the density of e-h pairs. 

The phenomenological equation ^ for study of periodical structures in exciton systems was 
employed in J2| for the case of spatially uniform exciton generation in bulk materials. In the 
case of infinite exciton lifetime (excitons are neither created nor decay), the above equation 
describes equilibrium states, i.e. two phases with the densities 

ntc = <x ± 1 (9) 

(in dimensionless units). The solution determines the density of the gas phase, while 
stands for the condensed one. In the non-equilibrium system with a finite value of the exciton 



Ordered Dissipative Structures in Exciton Systems in Semiconductor Quantum Wells 



5 




Figure 2. Exciton density nexix,y) for the parameters corresponding to the generation rate given in 
Fig. 1: a) for the curve 1, b) for the curve 3. 



lifetime Tex, the stationary uniform solutions of equation (jUJ are unstable. The exciton density 
varies in space periodically in the range limited by certain values rimin and nmax- In the case 
of the finite lifetime nmin and rimax can differ essentially from their equilibrium values. In the 
limit of the long exciton lifetimes (rex — *■ oo) nmin and nmax approach, respectively, their values 
in the gas and the condensed phases. 

In the given contribution the generation rate G is an essentially inhomogeneous function of 
coordinates. The system of the equations (j^J, (|HI) is a closed system of the equations, which 
enables describing phase states in the exciton system under inhomogeneous excitation. 

4 Numerical simulation of exciton condensed phase formation 
on the ring 

4.1 Exciton concentration distribution 

Having determined the exciton generation rate by solving equation (jlj, we numerically studied 
the stationary solutions of the nonlinear Eq. Q for the exciton density. We chose the following 
parameters of the system: D^x = lOcm^s"^, Tcx = 10"'^ s, ^ = 4//m, r = 3 • 10~^ s, n^^ = 
1.875 • 10^ cm-2, = 1.2n°x, T = 2K, /3 = 0.015. The simulation showed that for the exciton 
generation rate of equation (0J) in Fig. 1, the excitons condense into a ring which can break down 
into separate fragments (see Fig. 2a). There were 28 and 36 islands on the rings of the radii 
242// m and 295/.i m for the curves 1 and 2 of Fig. 1, respectively. An island shape is oval (of the 
order of 20// m x 30// m at the half- maximum) , which is determined by the parameter K. 

Strictly speaking, in the system with strongly non-uniform density distribution, it is impos- 
sible to assign certain phase states to separate areas, nevertheless points of high density can be 
attributed to the condensed phase, and points of low density to the gas phase. Thus, the created 
fragments are periodically positioned islands of the condensed phase of excitons, emergence of 
which is caused by the finite value of lifetime and the inter-exciton attraction. The condensed 
phase comprised by a periodic array of exciton islands corresponds to the fragments observed 
in experiments |21 13 El • 

The radius of the ring of the condensed phase grows in two cases: 

1) with increasing the pumping; 

2) if the laser spot expands leaving the integral intensity constant. 

With expanding the ring radius, the number of the exciton islands increases. They remain 
of the same size, and the distance between them is practically invariant. Similar trends such as 



6 



A. A. Chernyuk and V.I. Sugakov 



increased radius of the exciton ring with increasing the irradiation intensity and in the case of 
the expansion of the excitation spot were observed in |21 El S El El El • 

Some experiments |11 El report observations of an inner ring in the photoluminescence 
emission. We attribute emergence of the internal ring to specific character of the irradiation 
intensity profile in the vicinity of the laser spot edge. If the exciton generation near the laser 
spot does not change abruptly then, as our study shows, an internal concentric ring may appear 
around the central maximum. Our calculations show that it can be fragmented as well. Yet, 
we have to note that application of our theory is limited by assumption that the exciton gas 
is believed to be cooled, which might not be true near the laser spot. Possible reasons of the 
inner ring appearance due to the non-homogeneous exciton distribution and temperature are 
stated in the paper ^^I- Additional investigations at different pumpings are needed for the final 
clarification of the inner ring formation. 

As already specified, with increasing the e-h recombination rate W exciton generation rate 
becomes sharper (compare the curves 1 and 3 in Fig. 1). The results of the calculations show 
that with increasing W the fragmented ring of the condensed phase becomes continuous (see 
Fig. 2b). This is seen from the comparison of Figs. 2a and 2b, obtained for the different values 
of W (3cm^s~^ and 22cm^s~^, respectively) and for the same values of other parameters. 

With decreasing K, the island size of condensed excitons diminishes, so the fragmented ring 
passes to continuous one. Increase of K gives rise to continuous ring formation, because islands 
grow. The fragments do not form because the surface energy is smaller in the case of the 
cylindrical shape of a fragment. As K is related with the surface energy, the role of the surface 
energy with K changing is not defining and the solid circle forms in the system. 

The equation ^ can have solutions both as a continuous ring and as periodically placed 
islands of condensed excitons, and the choice of one of them can depend on the inhomogeneity 
of the system or on the boundary (initial) conditions. 

4.2 Localized spots 

The appearance of the "localized bright spots" can be attributed to non-uniformities of different 
nature in the QW Oj^. We have undertaken an attempt to simulate them in our approach, 
assuming that there is a certain region that differs from the surrounding regions by quality 
of the contact between the well and the matrix in which the wells are grown. We assumed 
that the value of Tg, which characterizes the time of the establishment of the electron density 
equilibrium between the wells and the matrix, varies in space and expressed this non-uniformity 
in the vicinity of the point tq = {L, 0) as 



We carried out calculations of the exciton density profile governed by equation (jSJ . An additional 
island located in the region of the non- uniformity emerges in this case (see Fig. 3). Obviously, 
in the case of several defects, several islands would appear as observed in O El. This is 
a consequence of non-linearity of the system. So, inhomogeneity, which arises in the crystal as 
a result of technological processes, can be a nucleus of new fragments. 

Similar calculations can be performed for inhomogeneities of a different nature. So, the 
presented theory can describe condensed states localized near defects. 

4.3 Modelling of two spots excitation 

We also carried out a simulation of the excitation of the system by two irradiation sources for 
description of the experiment with two spatially separated lasers (31. We solved the equations Q 




(10) 
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Figure 3. Emergence of a localized spot in the exciton density distribution at a point of local non- 
uniformity given by equation (|10|1 with parameters B = 150, L = 190/im, c = 12.2/im. 
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Figure 4. Exciton density distribution in the system irradiated by two laser spots. The distance between 
the spots centers are: a) 84// m, b) 61// m, c) 42// m, d) 21// m. The pumpings are smaller than in Fig. 2. 



and (jSJ with two laser sources with the radial shape (jHI each. Our calculations show non-linear 
contributions of the excitation from two spots: rings become extended in mutual directions and 
at a certain stage of spots approaching merge into a common oval-shaped ring (see Fig. 4). So, 
in such a way the rings can interact. 

The results in Fig. 4 confirm the non-additive character of the rings interaction, observed in 
the experiment [Hj. 

4.4 Behavior versus temperature 

We studied |3] behavior of the system versus temperature. We assumed that the parameter a 
(see equation 0) can be represented according to the Landau model: 

a = a{T,-T), 

where a < and Tc is the critical temperature. With temperature rising, the fragmentation of 
the ring disappears, its intensity falls, the radial structure of the exciton density becomes dim 
(see Fig. 5). Such transition with temperature growth is observed also in [2]. 
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Figure 5. The radial profile of the exciton density ncx(p) at different temperatures: 1) 2K, 2) 1.6 K, 
3) 2.7 K, 4) 4.4 K. = 3.7 K. The curves 1 and 2 represent the fragmented ring, the curves 3 and 4 are 
the profiles of low continuous rings. 

5 Conclusions 

Thus, the results of our study of the exciton condensation in a coupled QW around the laser 
spot can be summarized in the following conclusions: 

1. Fragmentation of the external ring of exciton density occurs at a certain threshold value 
of the exciton generation rate (pumping). 

2. Transition from a fragmented ring to a continuous one occurs at increased value of the 
electron- hole recombination rate (i.e. for a narrower barrier between the QWs), with re- 
duction of the surface energy (the parameter K) and at higher temperature. 

3. In a QW with a macroscopic defect a localized island of the condensed phase may emerge 
inside the external ring. 

4. Possibility of the development of an internal ring depends on the character of decrease of 
irradiation intensity at the edge of the laser spot. 

5. Exciton luminescence rings caused by two spatially separated laser spots attract and even- 
tually form a common oval-shaped ring. 

The considered periodical structures exist only due to the external excitation and are the 
development of self-organization in the system of exciton condensed phases. According to the 
classification of Jl] , structures of such type which exist exclusively in non-equilibrium conditions 
are called dissipative. 
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